CMB Polarization Systematics Due to Beam Asymmetry: 
Impact on Inflationary Science 
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Cosmic microwave background (CMB) polarization provides a unique window into cosmological 
inflation; the amplitude of the B-mode polarization from last scattering is uniquely sensitive to the 
energetics of inflation. However, numerous systematic effects arising from optical imperfections can 
contaminate the observed B-mode power spectrum. In particular, systematic effects due to the cou- 
pling of the underlying temperature and polarization fields with elliptical or otherwise asymmetric 
beams yield spurious systematic signals. This paper presents a non-perturbative analytic calculation 
of some of these signals. We show that results previously derived in real space can be generalized, 
formally, by including infinitely many higher-order corrections to the leading order effects. These 
corrections can be summed and represented as analytic functions when a fully Fourier-space ap- 
proach is adopted from the outset. The formalism and results presented in this paper were created 
to determine the susceptibility of CMB polarization probes of the primary gravitational wave signal 
but can be easily extended to the analysis of gravitational lensing of the CMB. 

PACS numbers: 98.70.Vc 



INTRODUCTION 



Upcoming cosmic microwave background (CMB) po- 
larization experiments are poised to immensely improve 
our understanding of the early universe. The significance 
of polarization lies not only in the fact that it increases 
the amount of data provided by temperature anisotropy 
alone but also because it is more sensitive to several phys- 
ical processes which took place in the early universe, 
e.g. the gravitational wave background produced dur- 
ing inflation (Seljak & Zaldarriaga [1], Kamionkowski, 
Kosowsky & Stebbins [2]) and reionization (Zaldarriaga 
[3], and Fan, Carilli & Keating [4]). Partially (but not 
only) due to this fact, CMB polarization has a unique 
dependence on the basic cosmological parameters. This 
feature can be used to remove some of the degeneracies 
afflicting cosmological parameter estimation from tem- 
perature anisotropy alone. One of the CMB polariza- 
tion's main features is the dependence of its B-mode 
(curl-like polarization) on the amplitude of the stochas- 
tic gravitational wave background generically predicted 
by inflationary models. Due to its faintness, the B-mode 
polarization is prone to degradation by various system- 
atic effects on a wide range of scales and it is important 
to remove these spurious contributions. This must be 
done in addition to controlling the various systematics 
induced by diffuse galactic emission (Amblard, Cooray 
& Kaplinghat [5]) and contamination from E-B mixing 
due to partial sky coverage (Lewis, Challinor & Turok 
[6], de Oliveira-Costa & Tegmark [7], Brown, Castro & 



Taylor [8]). 

This paper describes an analytic approach to assess 
systematics induced by imperfections of the polarime- 
ter's main beam. Several effects which contaminate the 
power spectra stem from temperature and polarization 
variations over scales comparable to the beamwidth (Hu, 
Hedman & Zaldarriaga [9], Ponthieu [10], Rosset et al. 
[11], O'Dea, Challinor & Johnson [12]); these effects can 
be modelled and characterized by the spurious C; they 
produce. 

Our results are power spectra presented as a combina- 
tion of the underlying power spectra with mixing coeffi- 
cients which are infinite sums of analytic functions. In 
practice however, this infinite series must be truncated 
and therefore our result is effectively equivalent to a series 
expansion. As expected, for small beam imperfections we 
find that the higher order corrections contribute negligi- 
bly to the spurious polarization. We define the small 
parameters characterizing the systematic effects in Table 
I. They are; the gain factor g, the differential beamwidth 
of the beams /x, the differential pointing p, the beam el- 
lipticity e and the beam rotation e. 

This paper is organized as follows; in section 2 
we present the basic mathematical formalism of spin- 
weighted fields used to characterize the systematic ef- 
fects, beam convolution, and our analytic results for the 
temperature-polarization leakage and polarization con- 
version in a single beam. Ultimately, we consider bolo- 
metric polarimetry which is most conveniently described 
using the Stokes parameters (e.g, Masi et al. [13]) /, Q 
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and U (V = for CMB polarization) as opposed to the 
Jones matrix formalism which is particularly useful for 
describing coherent polarimeters (Hu, Hcdman & Zaldar- 
riaga [9], O'Dea, Challinor & Johnson [12]). These pa- 
rameters are derived from differences in intensity in the 
Gaussian 2-D polarized beam response function for each 
polarization. Our two-beam experiment and the induced 
power spectra are derived in section 3 and the numeri- 
cal results for the B-mode power spectrum are described 
in section 4. We end with a discussion of the impact of 
these effects in Section 5. 



MATHEMATICAL FORMALISM 

We work entirely in Fourier space and begin with the 
expansion of both temperature and Q and U Stokes 
parameters in plane waves. Since all the effects con- 
sidered here are due to the beam asymmetry, and the 
beamwidths are typically on the degree or sub-degree 
scale, we can safely employ the fiat-sky approximation as 
far as effects related to the beam shape and size are con- 
sidered. The underlying physical power spectrum is cal- 
culated by CAMB (Lewis, Challinor & Lascnby [14]) us- 
ing the full sky. The temperature and other (integer-spin) 
combinations of the Stokes parameters are expanded in 
harmonic space as follows (e.g. Zaldarriaga & Seljak [15]) 

(Q±iU)(n) = ^2a± 2 ,im ±2Yi m (n) (1) 

l,m 

where the expansion coefficients of the spin ±2 polariza- 
tion parameters can be presented in terms of the E and 
B polarization modes 



a±2,lm = Ei m ± iBi m . 



(2) 



E and B are scalar and pseudoscalar under parity (having 
even and odd parities) respectively, and are sometimes 
referred to as the 'electric' and 'magnetic' (or 'gradient' 
and 'curl') polarization components. In the flat sky ap- 
proximation, Eq.(l) becomes 



T(x) = j- I T,e a x d 2 l 



X±(x) = (Q±iU)(x) 



„- / {E l ± 2B 1 )e 4l - x e ±24( *-^ ) d 2 l (3) 

2-7T J 



where T\, E\ and B\ are the Fourier components which 
are functions of the wave vector 1 only. 

Since in real space the temperature and polarization 
patterns are convolved with the beams, these expres- 
sions are simply the product of their Fourier transforms 



in Fourier space. We restrict the discussion to an ellip- 
tical gaussian beam (with major and minor axes a x and 

BM 1 cxp ( {x ~ Px)2 fo-^H (4) 

and its Fourier transform is 



£?(!)= cxp JL» + »1 . p . 



(5) 



The pointing error merely shifts the phase of the beam 
representation in Fourier-space. It is useful to switch to 
polar coordinates at this point 

l x = I cos(0i + ip — a) 

l y = I sin(0i + ip — a) 

p x = pcos9 

p y — ps'm8 (6) 

where the angles ip, a and 6 are defined below. The 
Fourier representation of the beam (Eq.5) then becomes 

B(\)d 2 \ = e -f-^os2(^+V-a)+^cos(0 ! -a-e+V) WW ^ (7) 

where 

I 2 



I 2 



(8) 



The definitions of the parameters in terms of the mean 
beamwidth a, differential beamwidth p, and ellipticity e, 
are given in Table II. Employing the expansion of 2-D 
plane waves in terms of cylindrical Bcssel functions 

n— 00 

e a P cos( l p l -4> P ) = ^ i n J n (lp)e m( ' l "~' l >>>\ (9) 

n— — 00 

the definition of modified Bcssel function 



I n (z) = i n J n (iz), 
and the symmetry relation 

J- n (z) = (-l) n J n (z), 

Eq.(7) becomes 



(10) 
(11) 



OO OO 



B(l) = E l 2m+n Irn(z)Jn(lp) 



7i— — 00 m— — 00 



^ e i(2m-\-n)ip—in9 e i(2m-\-n)((l)i — a) 



00 00 



= J2 J2 B m . n e^ 2m+n ^-^ (12) 



n— — 00 m— — 00 



where a = f3 + 8 + is the angle of the polarization axis 
in some coordinate system fixed to the sky (Fig. 1). We 
will employ this relation repeatedly in this work. 
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FIG. 1: The beam profile with the angles /3, 6 and ip shown for a single beam. The horizontal x and vertical y axes are fixed 
to the focal plane. The axis of polarization sensitivity makes an angle ip with the ellipse major axis. 
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TABLE I: Definitions of the parameters associated with the 
systematic effects. Subscripts 1 and 2 refer to the first and 
second beams, respectively. The differential gain parameter 
g and the beam rotation parameter e are not related to the 
beam shape and therefore are not defined in Eqs. (7)-(8), but 
rather are global parameters which define beam mismatch. 
We defer the exact definitions of these parameters to the rel- 
evant sections. 



where the sums run over all measurements of the pixel 
p, A is given by 



Aj = (1, \e 2 



' 2 



(15) 



and A T is A transposed. In general the matrix elements 
of A depend on the angle a = (3 + 9 + ip + e where the 
angles 9 and ip were defined above, the angle e is the sum 
of uncertainties in these two angles, and [3 is the angle 
between some arbitrary axis and the x axis of the focal 
plane. Clearly, the angular coverage uniformity of each 
pixel depends on the details of the scanning strategy. We 
address this issue in the following section. 

General Considerations 



DUAL POLARIZED BEAM EXPERIMENT 

We represent the data collected by the detectors as 
(e.g. Tegmark [16]) 



d(p) = A(p)m(p) + n(p) 



(13) 



where d is the measured data, A(p), m(p) and n(p) are, 
respectively, the pointing matrix, map vector with coef- 
ficients T, Q — ill, Q + ill, and the noise, at the pixel 
p. For gaussian white noise, the optimal map to assumes 
the form (Tegmark [16]) 



m(p)= \J2 A J A 

liep / Vj'ep 



(14) 



The effects we consider in this work arise either from 
circular beams with unmatched main-beam full width at 
half maximum (FWHM; called the monopole effect) or 
from beams with differential ellipticities (quadrupole ef- 
fect). The effect of beam ellipticity on the temperature 
anisotropy power spectrum was considered earlier by, e.g. 
Souradeep et al. [17], and the effect of beam asymme- 
try on the two-point correlation functions of the tem- 
perature anisotropy and polarization was considered by 
Ng [18]. Our present work, however, is multipole space 
based. Also, differential gain or differential pointing 
(dipole effect) can induce further spurious polarization 
signals from temperature leakage due to beam mismatch 
as will be discussed below. Finally, if the polarization 
sensitivity axes of the beams are rotated (differential ro- 
tation effect) we expect mixing between the polarization 
E and B modes, and associated leakage of power between 
E and B. 
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We first consider a general case with e = (i.e., no 
rotation error). In this case, Eq.(14) for the two beam 
experiment reads 

I'(p) - [{A T A) + R,{A T A)R-% 1 

x [(A T A)I 1 (p) + R n (A T A)R- 1 I 2 (p)] (16) 

where the vector I(p) = [T,Q — ill, Q + iU](p), the sub- 
scripts 1 and 2 refer to the first and second beams, re- 
spectively, and i?,r is the rotation matrix by w. Using 
Eq.(15), we obtain 



T' = (T + ) + -((Q_-z[/_)e 2 -) 
+ \{{Q-+iUJ)e-*°) 
Q'±iU' = l(Q + ± + A (T_e ±2M ) 



(17) 



D 
1 



D 



-<(Q +T £/ + ) e ± 4 ») 



1 

D 
2 



{Q + TiU + ){e^ a ) 



where 



- j-(T_e T2la )(e ±4la ) 



X ± = -(X 1 ±X 2 ) 



(18) 



(19) 



with X = T, Q or U and D = 1 - (e 4ia ) (e" 4ia ) . Angular 
brackets stand for averaging the trigonometric functions 
due to the scanning strategy To simplify the following 
discussion we encapsulate the properties of the scanning 
strategy 

f(m,n) = ( e -i(2m+«)«) 
h±(m,n) = l[/(m, n) - f(m± 2, n){e ±4ia )]. (20) 

As an example of / and /i±, for an ideal scan strategy, 
the angle a is uniformly sampled, and as can be seen 
from Eq.(20), only combinations of n and to which satisfy 
2m+n = result in non- vanishing /(to, n) and h±{m, n). 
When the scanning strategy is non-ideal, and the beam 
mismatch does not have the required quadrupole symme- 
try, the above functions f(m,n) and h±(m, n) couple to 
the beam mismatch to satisfy the necessary properties of 
spin ±2 fields for certain combinations of to and n (see 
Eq. 27 below). Employing Eq.(3) we obtain the tem- 
perature anisotropy and polarization modes in Fourier 



space 

m.n 

+ \ Wi,« e * (2m+ " )0; (^1 - ^0] * Mm, n) 

m,n 

+ \ J2[( B -)m-i, n e i{2m+n)4 " {Ei + iBi)] * Mm, n) 

m,n 

E[±iB[ = Y,[{B+) m , n e K2m+n) * l (Ei ±iBi)] *ti ±l (m,n) 
+ J2i( B +)m±2^ 2m+n) ^ (Eh T iBi)} * h±i(m, n) 

m.n 

+2Y J [{B-) m ±i, n e l{2m+n) ^fi] * h ±l (m, n) (21) 



where * stands for convolution in Fourier space and 
(-B)m,n are defined in Eq.(12). 

The systematic power spectra are obtained by taking 
the angular average (in Fourier space) of the squared 
modulus of the expressions in Eq. (21). Due to the as- 
sumed statistical isotropy of the underlying power spec- 
tra, the celestial temperature and polarization can be 
taken out of the integral 



c[ XY = (xy*) = /xr? 

J 2ir 



(22) 



and we are left with integrations over the beam profiles 
in Fourier space. We use the auxiliary quantities 

A ={Ti(E-iB)*} 
C± = {(E + iB)(E ± iB)*) (23) 

in terms of which the power spectra are 





= iRc(C+ + C_) 


C[ TE 


= Rc(A) 


C[ B 


= iRc(C + -C_) 


C[ TB 


= lm(A) 


(J'EB 





and 



C[ L = {\T'i\ 2 }. 



(24) 



(25) 



The explicit forms of A and C± are given in the Ap- 
pendix. 

Two special cases, for which these results simplify con- 
siderably, are an ideal isotropic scan (in which every pixel 
is being scanned many times in a random orientation so 
both /(to, n) and h±(m, n) identically vanish (except for 
the case 2m + n = 0) as do their Fourier transforms /(l) 
and h±(l)), and uniform coverage of the observed field 
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(the number of hits per pixel is not large so / and h± 
get constant, non-zero values). In this case the Fourier 
transforms of / and h± are <5-functions and there is no 
multipolc- mixing in Eq.(21). Clearly, the first case is the 
limit of the second when A^hits is very large. In principle, 
the field of view can be few angular degrees so one may 
wonder if the convolution in Eq.(21) in Fourier rather 
than harmonic space is warranted. Indeed, all the un- 
derlying power spectra peak at high multipoles and the 
effects we consider here are on scales of few percent of 
the beamwidth which is assumed to be 1° at most and 
if the sky is covered relatively uniformly, the supports of 
the functions /(l) and h±(l) are narrow with a support 
on a very small range of Al. Combining these two facts 
together it is evident that we can still work in the flat-sky 
approximation in most cases of interest. 

We limit the following discussion to gaussian beams 
but similar calculations can be readily done in the case 
of other beam shapes (at least numerically). The cor- 



parameter 


beam 1 


beam 2 




<r(l + /x)(l + e) 


a(l-p)(l + e) 


ffy 


ff(l + /i)(l-e) 


«r(l-/i)(l-e) 


y 


^(l + M ) 2 (l + e 2 ) 


(Ml { l-p) 2 (l + e 2 ) 


z 


(la) 2 (l + p) 2 e 


(la) 2 (l-p) 2 e 



TABLE II: Definitions of the parameters associated with the 
dual-beam experiment. 



rections to the underlying power spectra are defined as 
follows 

ACf = C' l z (a,g i ,iH,e i ,p i ,e)-Ci z (a,0) 

Z e {TT, TE, EE, BB, TB, EB} (26) 

which are functions of the small parameters g, p, p, e, 
e and the underlying power spectra, where here Z = 
TT, TE, EE, BB, TB and EB, and as we see from Eqs. 
(21) and (24), the new power spectra are, in principle, 
combinations of all the underlying power spectra. Also, 
it is worth mentioning here the coupling between the var- 
ious systematic effects as far as higher order corrections 
are concerned. While Eqs. (21) and (12) are presented in 
terms of formally-exact, well-known functions, for calcu- 
lations the infinite series (Eq. 12) must be truncated, and 
depending on the degree of asymmetry and the physical 
scale in question (tantamount to the multipole number I), 
the number of terms in the scries considered determines 
the accuracy of the calculation. The scalings of the lead- 
ing order terms of the effects considered here are given 
in Tables III- VI assuming the underlying sky is unpolar- 
ized (except for the effect of rotation, to be discussed in 
the next section) . In obtaining these expressions we also 
assumed the scanning strategy is statistically isotropic, 
a reasonable assumption that significantly simplify the 
scaling relations we obtain. We have used the following 



definitions 

h = ^IM-i,o)| 2 

h = i|M-i,-i)l 2 + ^IM-M)l 2 

h = \(mi)h*_ (i,-i)) (27) 

where the functions f(m,n) and h±(m,n) are defined 
in Eq.(20). Since the leading orders of both J n (z) and 
I n (z) arc oc z n it is clear that terminating the series at 
some order n is equivalent to a power series expansion 
exact up to this order and not higher. Therefore, the 
calculation described here actually amounts to a power 
series expansion in the most general case, but simplifies 
considerably when cither the dipole or quadrupole effects 
can be ignored. Tables V and VI show the scaling rela- 
tions in the case of ideal scanning strategy corresponding 
to Tables III and IV, respectively. These are the irre- 
ducible signals that persist irrespective of the scanning 
strategy. One further remark is in order here: the cou- 
pling between the various effects is important when the 
parameters l 2 a 2 p, Ip and l 2 a 2 e are not negligibly small 
compared to 1; in this case higher order corrections are 
required and the cross terms that include the coupling 
between the effects cannot be ignored, but this is seldom 
the case. For given parameters p, e and p, there is a 
'critical' value of I beyond which higher order terms be- 
come important and in this range our calculation may be 
particularly useful. For a « 1° and e, p and pja « 0.1 
higher order effects become important at I « 1000 which 
is tantalizingly close to the scale at which the B-mode 
signal from lensing peaks (Zaldarriaga & Seljak [19], Hu 
[20]). However, it happens deep beyond the beam di- 
lution scale, I « 200. Only when p is not very small 
compared to a or e and p arc not small compared to 
unity do these higher order corrections contribute to the 
spurious polarization. 

Rotations 

We now discuss the rotation error ignored in the above 
treatment. Overall rotation of the two beams (cross- 
polarization) mixes E and B-modes and even induces T- 
B correlation. Non-orthogonality of the beams can be 
described in two different ways. If only one of the two 
detectors is miss-oriented by 5, there is induced E-B mix- 
ing. This is the approach taken by Ponthieu [10]. If the 
two detectors are disoriented, one by +5/2, the other one 
by —5/2, the Q- U mixing induced by the first detector 
is compensated by that of the second detector and there 
is no E-B mixing. In the following, we focus on the lat- 
ter when we refer to non- orthogonality of the beams. We 
here derive the expected signal. For (overall) rotation 
error e and non-orthogonality measured by 5 (which we 
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split between the two beams; one is missoriented by (5/2 
and the other by —5/2) 



5> T d,)(p) =|| M || I 



(28) 



where I = (T, Q - iU, Q + iU) and 
II M || = 



/ 



1 



2ia\ 2ie+iS 1, 



Ue 2la )e 



-2ia 



-2is—i8 



l„-2te- 



1 „2ie+i<5 



(29) 



-2ie—iS 



and upon using Eq.(16), we obtain 




T+%sm5^2 ± F±{QTiU)e 
(Q - iU) cos 8e 2l£ 
{Q + iU) cos(5e- 24£ 

where F± = =f1. Employing Eqs.(3) 

Qi ± if/i = (Si ± iB x )e^ 2i ' t ' 1 

we obtain to leading order 



(30) 



(31) 



n 



Ei - 2eB 1 
Bi + 2eE 1 



(32) 



and therefore 



AC, 
AC, 
AC, 
AC, 
AC, 
AC, 



T ~ 

t' e' 

E' ^ 
I ~ 

B' ^ 



Ae 2 C? 



T'B' 



E B 



TE 



2sC, 
2eCp 



(33) 



i.e. the leakage is a small fraction of Cf but can still be 
a significant B-mode contaminant. 



Monitoring the Contamination 

A potentially useful diagnostic is the T — B cross cor- 
relation. As mentioned above, this correlation function 
vanishes in the standard model. The spurious effects dis- 
cussed in this paper cause T to leak to both E and B, and 
therefore the correlations Cj B and Cf B do not generally 
vanish (see Eqs. 24 and Table IV). Furthermore, as we 
will see below, beam rotation can also induce Cj B and 



Cf B because rotating the telescope by e is indistinguish- 
able from rotation of the polarization plane by —e. As a 
result, power leaks from T-E to T-B and from E-E to E-B. 
A nonvanishing Cj B may be attributed to an imperfect 
removal of the spurious polarization signals discussed in 
this paper. However, beam systematics are not the ex- 
clusive generating mechanisms of Cj B and Cf B . A few 
physical sources of parity-violating correlations were al- 
ready discussed in the literature; parity-violating terms 
in the lagrangians of the electromagnetic and gravita- 
tional sectors (e.g., Lue, Wang & Kamionkowski [21], Liu, 
Lee & Ng [22], Feng et al. [23], Saito, Ichiki & Taruya 
[24], Xia ct al. [25], Komatsu et al. [26]), Faraday ro- 
tation at last scattering (Kosowsky & Loeb [27]), and 
hypothetical primordial helical magnetic fields (Caprini, 
Durrer & Kahniashvili [28]). The systematic T-B and E- 
B correlations may interfere with this exotic physics and 
a careful analysis of these correlations in required in the 
presence of beam systematics (Shimon & Keating [29]). 



NUMERICAL RESULTS 

The Stokes Q parameters associated with the polar- 
ization due to the optical imperfections discussed in this 
paper are shown in Figure 2. We have calculated all pos- 
sible six power spectra (Eqs. 24 and 33) for two different 
average beam widths FWHM 56' and 5'. Figures 3-8 de- 
pict beam rotation, differential pointing and differential 
cllipticities effects on the B-mode power spectrum for 
average beamwidths of 56' and 5', respectively. The dot- 
dashed lines refer to the inflation-induced i?-mode from 
primordial gravitational waves with tensor to scalar ra- 
tios of T/S = 10" 1 , 10^2, 1(T 3 and 1CT 4 , respectively, 
where we have used the definition used in CAMB for the 
tensor-to-scalar ratio. For the plots of the second order 
differential pointing effect the specified pointing error p 
refers to one of the beams (we left the angle 9 of this beam 
a free parameter), while the other beam has been as- 
sumed to have no pointing error . The quadrupole effect 
was calculated assuming the two beams have the same 
specified cllipticities, |e|, and here we left the angles that 
the polarization axes make with the major axes of the two 
beams, ipi an d V'2, as free parameters. We have de- 
fined = cos(2V>i) — cos(2f/>2), s,/, = sin(2?/>i) — sin(2i/>2) 
and eg = cos 26, sg = sin 29. Note that in the case consid- 
ered here, that both beams have the same ellipticity e|, 
we obtain, as expected, that there in no spurious polar- 
ization when -01 = ip2- Similarly, had we assumed both 
beams have the same pointing p we would have obtained 
no spurious polarization if 9\ = 9 2 in the case of the ideal 
scanning, where second order pointing is the leading or- 
der contribution (not so in the case of non-ideal scanning, 
where the leading order is the dipole effect). 

As expected, the larger the beamwidths the larger the 
angular scales on which the systematics peak. The sys- 



FIG. 2: An illustration of the monopole, dipole, quadmpole and gain effects; Q parameter only is depicted. 



B-Mode Polarization (FWHM=56 arcmin), rotation 

10 2 | . . 1 1 1 . r- 




io- 8 - 



10° 10 1 10 z 10 3 

I 

FIG. 3: The contribution of differential rotation to the B-mode power spectrum (56' average beamwidth). Shown are the 
effects for e = 0.01, 0.02, 0.03, 0.04 and 0.05 of a radian. For comparison, the dot-dashed curves refer to the contribution 
from primordial gravitational waves with tensor to scalar ratios T/S = 10 _1 , 10 -2 , 10 -3 and 10 -4 . The dashed curve is the 
B-mode polarization produced by gravitational lensing by the large scale structure. 



tematics are calculated for the small parameters; g, /x and 
e = 1, 2, 3, 4 and 5 % of the mean beam and p and s the 
same fractions of the mean beam and a radian, respec- 
tively. The scaling relations of these effects are given in 
Tables III- IV (general scanning strategy) and V-VI (ideal 
scanning strategy). As explained in the last section, the 
higher order corrections (beyond the leading order) are 
important on some scale, typically I rj er /i — ' which 
for a » 1° and e and /i on the 1% level corresponds to 
I w 2000. On this scale, the S-mode from the CMB 
lensed by the large scale structure (e.g. Zaldarriaga & 
Seljak [19], Hu [20]) is non- negligible (the dashed lines 
in Figures 3 through 8). However a beam that size is 
insensitive to features on scales of I ~ 2000. It is also im- 



portant to mention here that the effect due to rotation 
(Figures 3 and 6) closely follows the C, shape and merely 
reflects the leakage of E to B due to rotation (Eq. 33). 
One more remark should be made regarding the coupling 
of the differential pointing and quadrupole effects men- 
tioned above. As seen from Eqs. (12) the pointing error 
and ellipticity are coupled through products of the cylin- 
drical and modified Bessel functions. In practice however 
the coupling is small since the pointing error parameter 
Ip is small. In this case, since the leading term of J n (x) 
is oc x n it is safe to consider only the leading terms in 
the infinite sum in Eq. (12), this approximation consid- 
erably simplifies numerical calculations. In calculating 
the plots we assumed perfect scanning strategy so the 



8 



B-Mode Polarization (FWHM=56 arcmin), pointing 
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FIG. 4: The contribution of differential pointing to the B-mode power spectrum (56' average beamwidth). The values shown 
should be multiplied by Sg (Table III). In case 8 = 0° all the spurious polarization is in the E-mode. For a given p, the maximum 
B-mode is obtained at 8 = 45°. Shown are the effects for p — 0.01, 0.02, 0.03, 0.04 and 0.05. For comparison, the dot-dashed 
curves refer to the contribution from primordial gravitational waves with tensor to scalar ratios T/S = 10 _1 , 10 -2 , 10~ 3 and 
10" 4 . The dashed curve is the B-mode polarization produced by gravitational lensing by the large scale structure. 



B-Mode Polarization (FWHM=56 arcmin), quadrupole 
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FIG. 5: The contribution of the differential beam ellipticity ('quadrupole' effect) to the B-mode power spectrum (56' average 
beamwidth). The values shown should be multiplied by s 2 (see Table III). Shown are the effects for e = 0.01, 0.02, 0.03, 0.04 
and 0.05. For comparison, the dot-dashed curves refer to the contribution from primordial gravitational waves with tensor 
to scalar ratios T/S = 10" 1 , 10" 2 , 10" 3 and 10" 4 . The dashed curve is the B-mode polarization produced by gravitational 
lensing by the large scale structure. 
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B-Mode Polarization (FWHM=5 arcmin), rotation 
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FIG. 6: The contribution of differential rotation to the B-mode power spectrum (5' average beamwidth). Shown are the effects 
for e = 0.01, 0.02, 0.03, 0.04 and 0.05 of a radian. For comparison, the dot-dashed curves refer to the contribution from 
primordial gravitational waves with tensor to scalar ratios T/S = 10 _1 , 10 -2 , 10 -3 and 10 -4 . The dashed curve is the B-mode 
polarization produced by gravitational lensing by the large scale structure. 



monopole and gain identically vanish as expected and as 
can be easily verified from Eqs. (12), (21) and the re- 
lations in the Appendix. Again, it is important to note 
that the power spectra for the second order pointing er- 
ror (Figures 4 and 7) and quadrupole (Figures 5 and 8) 
should be multiplied by the functions eg and sg, and 
and s^, respectively, as described in Tables III, IV, V 
& VI. These functions can vanish. For instance, the ef- 
fect of ellipticity will not contribute to the B-mode power 
spectra if the polarization sensitive axes are parallel to 
one of the ellipse principal axes (ip — 0,7r/2), there is no 
shear of the held in this highly symmetric case and all 
the spurious polarization is in the E-mode as can be seen 
from Table V. This is no longer the case with non-ideal 
scanning strategy (Table III); in this case the hrst order 
dipole do not depend on 9. Also, if ipx = -02 there will 
be no induced polarization by differential ellipticity, not 
in the E-, nor in the B-mode. Similarly, for the pointing 
error, if 9 = 0, all the spurious polarization will con- 
taminate the E-mode (if the scanning strategy is ideal). 
The power spectrum associated with the underlying sky 
was calculated by CAMB using cosmological parameters 
consistent with WMAP (Spergel et al. [30]). We ignore 
gravitational lensing and the tensor contribution to the 
underlying sky, but we do show them in the figures for 
reference. We stress again that nonvanishing C[ TB and 
C[ EB may indicate an incomplete removal of the spurious 
polarization signals discussed here. The relevant expres- 
sions are given in Tables IV & VI. Figures 9 to 11 show 



that C[ TB peaks at the few /iK level (all these values are 
for 'maximum shear', i.e. 9 = 45° or i/j — 45°), and so 
constitute only upper limits for the given parameters p 
and e. Only the contributions from (second order) dif- 
ferential pointing, differential ellipticity and differential 
rotation for the TB cross-correlation are displayed since 
the monopole and gain contributions identically vanish 
for an ideal scanning strategy. One more aspect of this 
analysis is the higher order corrections of the effects stud- 
ied in this paper. We find that these corrections have lit- 
tle effect even on the smallest scales (where the gradients 
are large) in the ideal cases we studied. These effects 
become relatively large on scales much smaller than the 
mean beam size where beam dilution is significant. Wc 
conclude that at least for nearly ideal scanning strategy 
these higher order corrections can be safely neglected. 



DISCUSSION 

The systematics discussed here, more than simply af- 
fecting the peak of the B-mode power spectrum at I « 
100 are likely to impact the polarization signal due to 
lensing at I rj few hundreds as shown in Figures 3 to 
8, at least for small beams. A full analysis of the ef- 
fect of gravitational waves on CMB B-mode polarization 
is underway, and requires an assessment of the effect of 
gravitational lensing. The results presented in this pa- 
per show that the beam systematics are likely to further 
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B-Mode Polarization (FWHM=5 arcmin), pointing 
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FIG. 7: The contribution of differential pointing to the B-mode power spectrum (5' average beamwidth). The values shown 
should be multiplied by s|. Shown are the effects for p = 0.01, 0.02, 0.03, 0.04 and 0.05. For comparison, the dot-dashed 
curves refer to the contribution from primordial gravitational waves with tensor to scalar ratios T/S = 10 _1 , 10 -2 , 10 -3 and 
10" 4 . The dashed curve is the B-mode polarization produced by gravitational lensing by the large scale structure. 



B-Mode Polarization (FWHM=5 arcmin), quadrupole 
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FIG. 8: The contribution of the differential beam ellipticity ('quadrupole' effect) to the B-mode power spectrum (5' average 
beamwidth). The values shown should be multiplied by s^. Shown are the effects for e = 0.01, 0.02, 0.03, 0.04 and 0.05 
of a radian. For comparison, the dot-dashed curves refer to the contribution from primordial gravitational waves with tensor 
to scalar ratios T/S = KT 1 , 10~ 2 , 10~ 3 and 10" 4 . The dashed curve is the B-mode polarization produced by gravitational 
lensing by the large scale structure. 
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T-B Cross Correlation (FWHM=56 arcmin), rotation 



O 




FIG. 9: The contribution of differential rotation to the T — B cross correlation (56' average beamwidth). Shown are the signals 
for e = 0.01, 0.02, 0.03, 0.04 and 0.05 of a radian. Dashed lines correspond to negative values of the cross- correlation induced 
by the underlying T-E cross correlation. 



T-B Cross Correlation (FWHM=56 arcmin), pointing 




FIG. 10: The contribution of differential pointing to the T — B cross correlation (56' average beamwidth). The values shown 
should be multiplied by se- Shown are the effects for p = 0.01, 0.02, 0.03, 0.04 and 0.05. 



complicate this process both directly (with a residual ef- 
fect on degree scales) and indirectly (by mimicking the 
effect of lensing). While the lensing signal can be par- 
tially extracted by invoking optimal estimators that use 
the non-gaussianity of the lensed anisotropy and polariza- 
tion (Hu & Okamoto [31]), the systematics discussed here 



cannot be removed by this method because they inherit 
the statistics of the underlying sky, which is gaussian in 
the standard model. The lensing-induced polarization is 
a direct probe of structure formation processes and is 
sensitive to few cosmological parameters, most notably 
the neutrino mass. An accurate analysis should take the 
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beam effects discussed here into account. Our calcula- 
tion demonstrates that at most, these systematics are on 
the level of few tens to few hundredths of nK and their 
significance as contaminants of the primordial B-mode 
depends on the tensor to scalar ratio. Comparing our 
results to the results of Ponthieu [10] it is evident that 
the spurious power spectra calculated here are less noisy 
because they do not include the cosmic variance. The 
formalism described here directly employs the underly- 
ing power spectrum as opposed to Ponthieu [10] whose 
results are based on producing synthetic CMB maps and 
their first and second spatial derivatives with HEALPix 
(Gorski et al. [32]). These are combined with, e.g. the 
first order moments of the effective beam and hits dis- 
tribution in each pixel, producing simulated 'observed' 
maps (including their systematic effects), and these maps 
can be studied by themselves or used to produce power 
spectra. Although our analytic results are given in their 
most general form, including the effect of scanning strat- 
egy, in practice /i, / 2 and f 3 (Eq. 27) should be numer- 
ically calculated. This is consistently done by the ap- 
proach employed by Ponthieu [10]. Our analysis shows 
that the contributions from higher order corrections and 
coupling between the various effects may be of some im- 
portance in principle, but for the cases studied here of 
a perfect elliptical gaussian beam with mean beamwidth 
rs 1° and cllipticity and pointing error on the 5% level 
with a perfect scanning strategy, we find that the higher 
order corrections have a negligible effect on the angular 
scales of interest. Another aspect is the potential use of 
the nonvanishing Cf B and Cf B power spectra as moni- 



tors of an imperfect removal of the spurious beam effects 
using the fact that these two power spectra are expected 
to identically vanish in the standard cosmological model. 
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APPENDIX A: AUXILIARY CORRELATION 
FUNCTIONS 

We list here the auxiliary correlation functions A and 
C± 

C± = {(E + iB)(E ± iB)*) = 

(h+{m, n)h*±{m', n ') e <P("»-"»')+(»-«')]*«) 

mnm'n' 

* [^CT{B-)m+l,n{B-)* m ,± 1>n , 

+ 2Cf E {(B-) m+ i^ n ((B+) m > n i + {B + ) rn i±2,n'T 
+ [B-)m'±l,n> (( B +)m,n + {B+) m+ 2,n)] 



+ 


Q [{B+)m'n' + { B +)*m'±2,n>] 




X 


[{B+)m,n + (B+) m +2,n] 




± 


Q S ((B+)m,n — (-S+)rn+2,n) 




X 


((B+)m>,n> - (B+)m>±2,n>) 




± 


2iC ; TB [(i?_)„ l+ i ! „ ((B + ) m i±2,n' 


- (-B+)m' 


± 


( B -)*n>±l,n> ((B+)m,n - { B +)-m- 


t-2,ra)] 


4- 


■/~lEB\/ 1 T3 \ , (r> \ \ 
)m,n + {-&+ Jm+2,n) 




X 


(\B+) m >±2 n' ~ v®+)m',n') 




± 


\{B+)m,n ~ {B+),n+2,n) 




X 


(( B +)m'±2,n' + ( B +)* m >,n>)] 






/TV 771 ■ r>\*\ 

{1 {h ~iB) ) = 






2 y {J(m,n)h_(m ,n )e 


m / )+(n— ri 








+ 


Q T [(^+)mn ((-B+)m'n' + 


<7i' — 2,n' ^ 


+ 


( 5 -)m'-l,n' ((■S_) m +l,n + (#-) 


m — 1 , n ) ] 


+ 


2^/ B [(-®+)m'n' + (^+)m'-2,n'] 




X 


[(-B_) TO +l )T j + (i?-) m -l,n] 




+ 


2 C 'i B [(- B +)m',n' - ( B +)*n> -2,n>] 




X 


\{B-) m +l,n — {B-) m -l,n\ 




+ 


iCf B [(B + ) mn ((B + ) m i <n i — (B + 


)m'-2,n') 


+ 


{ B -)*n'-\,n' ((B-)m-l,n - (B-) 


m+l,n)j 


+ 


—C EB [((B-) m+ i^ n + (B_) m _i )? 





X 


((B+) m >n> - (B+)*„i-2,n>) 






((-B-)ro+l,n — (-B-)m-l,n) 




X 


((- B +)m'.n' + (' B +)m'-2,n')]] 





(A.l) 



(A.2) 

where * denotes 2-D convolution. 
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effect 
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<? 2 /i*CT 


monopole 
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4/i 2 (/a) 4 /i*Q T 


4^ 2 (Za) 4 /i*Q r 


pointing 
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±c»Cf [1 + J (lp)]J 2 (lp) 




slCjJl{lp) 






-ceJ?{lp)C?*f 3 


+j 1 2 (; P )c i T */ 2 


-Jl(lp)Cj*h 


quadrupole 
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-h(z)h{z)c^CT 




ll{z)s%CT 


rotation 







^c B 


4 e 2 Cf 



TABLE III: The leading order contributions of the systematic 
effects to the power spectra Cf B , Cf and Cf assuming the 
underlying sky is not polarized (except for the rotation signal 
where we assume the E-, and B-mode signals are present) and 
general sky scanning. 
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2eCj E 
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TABLE IV: The contribution of the systematic effects to the 
power spectra Cf B , Cf B assuming the underlying sky is not 
polarized (except for the rotation signal when we assume E-, 
and B-mode polarization are present) and general sky scan- 
ning. 
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TABLE V: The contribution of the systematic effects to the 
power spectra Cj B , Cf and Cf assuming the underlying 
sky is not polarized (except for the rotation signal where we 
assume the E-, and B-mode signals are present) and ideal sky 
scanning. 



effect 


parameter 


ACf B 


ACf B 


gain 


9 








monopole 


A* 








pointing 


P 


±[1 + Mlp)]J2{lp)SgCT 


SQCgJl(lp)CT 


quadrupole 


e 


-i a { z )h{z). H cT 




rotation 


£ 


2£CJ E 


2eCf 



TABLE VI: The contribution of the systematic effects to the 
power spectra Cf B , Cf B assuming the underlying sky is not 
polarized (except for the rotation signal when we assume E-, 
and B-mode polarization are present) and ideal sky scanning. 



